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In this paper we establish a formalism for the computation of observables due to acceleration- 
induced particle physics processes. General expressions for the transition rate, multiplicity, power, 
spectra, and displacement law of particles undergoing time-dependent acceleration and transitioning 
into a final state of arbitrary particle number are obtained. The transition rate, power, and spectra 
are characterized by unique polynomials of multiplicity and thermal distributions of both bosonic 
and fermionic statistics. The acceleration-dependent multiplicities are computed in terms of the 
branching fractions of the associated inertial processes. The displacement law of the spectra predicts 
that the energy of the emitted particles is directly proportional to the accelerated temperature. 
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I. INTRODUCTION 

As a first order approximation to quantum gravity, quantum field theory in curved spacetime predicts the creation 
of particles from the vacuum. Of particular importance is the creation of thermalized particles via the Parker [1], 
Hawking [2], and Unruh [3] effects. These particle production mechanisms have given considerable insight into the 
quantum properties of spacetimes which contain horizons. The Rindler horizon present in accelerated reference frames 
enables one to use the Unruh effect, and its associated thermal bath, to induce particle transitions, i.e. decays and 
excitations. Our current understanding of these processes begins with the original computation of Mueller [4] which 
showed that one could use acceleration to change the decay rate of weakly interacting particles as well as excite 
stable particles into more massive states. Matsas and Vanzella [5-7] extended the analysis to fermions, analytically 
computed the power emitted during the transition, and numerically computed the spectra of the final state particles. 
In a previous work [8], we used scalar fields to extend the formalism to encompass transitions with arbitrary final state 
multiplicity. In this paper, the analytic computation of both the power and the spectra for arbitrary multiplicity is 
carried out. We also import certain relevant inertial quantities, such as the branching fractions of various decays, into 
the formalism and compute the acceleration scale to select the relevant decay pathway, i.e. multiplicity. Moreover, 
with the Planckian spectra obtained, we compute the peak energy of the emitted particles via a generalization 
of Wien’s displacement law. This establishes that the most probable energy of the emitted particles is peaked 
about the accelerated temperature. These results are punctuated by the entire analysis being carried out using a 
newly developed time-dependent formalism which agrees with the previous developments of Obadia and Milgrom [9], 
Kothawala and Padmanabhan [10], and Barbado and Visser [II]. The time dependence and ability to compute a wide 
class of observables developed in this paper establish a basic foundation for an acceleration-induced particle physics 
phenomenology with applications to highly accelerated systems e.g. the inflationary epoch of the early universe or 
perhaps the evaporation of black holes. 

In this paper. Sec. II computes the generalized response function. Using a fully relativistic treatment, i.e. quantized 
fields, the response function is derived in a similar manner as one would in the case of constant acceleration. However, 
to extend the analysis to encompass time dependence, the generalized response function is expressed as the Fourier 
transform with respect to the rapidity rather than the proper time. This prescription allows for the use of the 
same mathematical treatment in computing observables as in the case of constant proper acceleration. Section HI 
focuses on the time-dependent formalism utilized to compute various spacetime quantities for use in characterizing the 
accelerated motion, e.g. trajectories, spacetime intervals, and Lorentz gamma. These quantities are used to evaluate 
the Wightman functions and its variants which correspond to the relevant observables. In Sec. IV, the appropriate 
variants of the Wightman functions corresponding to the transition rate, power, and spectra are computed for arbitrary 
trajectories. Each is then evaluated along the appropriate time-dependent trajectory derived in the previous section. 
In Sec. V we evaluate the acceleration-induced transition rate. In doing so, we also derive the transition rate 
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polynomials of multiplicity, characterized by an integer index, along with the associated bosonic thermal distribution. 
Section VI utilizes the transition rates to compute the acceleration scale that selects the dominant multiplicity based 
on the branching fractions of the relevant processes in the inertial limit. Specific limiting cases are considered which 
yield analytic results. The radiated power is computed in Sec. VII using the appropriate variant of the Wightman 
function. In doing so, we also obtain the power polynomials of multiplicity, characterized by a half-integer index, 
along with the associated fermionic thermal distribution. Section VIII contains the derivation of the energy spectra. 
The resulting Planck-like spectra are expressed in terms of the spectra polynomials of multiplicity, characterized by an 
integer index, and the associated bosonic thermal distribution at finite chemical potential. In Sec. IX the displacement 
law of the peak energy emitted is computed for arbitrary acceleration, multiplicity, and transition energy. The result 
generalizes Wien’s displacement law, in that the peak energy is proportional to the accelerated temperature weighted 
by a numerically determined parameter. Final comments and a summary of the conclusions are presented in Sec. X. 
All calculations are performed using the natural units h = c = ks = i- 

II. GENERALIZED RESPONSE FUNCTION 

In this section we set up a formalism capable of computing a wide class of particle transitions regardless of the 
number of final state products [4,8]. To facilitate the analysis, we consider all particles to be scalars. Consider an 
initial Rindler particle moving along an arbitrary time-dependent accelerated trajectory and decaying into a hnal 
state containing ur Rindler particles and n Minkowksi particles. This process is schematically written as 


'fl + 'i’2 + ■ ■ ■ + + <l>2 + ■ ■ ■ + (l>u- ( 1 ) 

The initial and final Rindler particles are denoted by and are used to describe any particle under acceleration. 
We consider the initial Rindler particle to be massive while the final state Rindler particles can be massless or massive 
in any combination. The massless Minkowski particles in the final state are denoted by (j)k and are used to describe 
any particles propagating along inertial trajectories. To describe these transitions, we consider their coupling to be 
described by the following general interaction action, 

^ p 'kl 

Si = ^ n n (2) 

^^ d r—l m—1 

The coupling constant G„ is labeled by the multiplicity and can be determined by fixing it to a known process 
in the inertial limit. We note that the coupling may be dimensionful depending on the transition in question. The 
term is included to absorb the overall normalization of our Rindler states. The domain of integration is confined 
to the right Rindler wedge where the modes of our accelerated fields are defined. We denote the Fock states of our 
Rindler particles [4'^) with the index j characterizing their energies. As usual, we label our Minkowski states |k^) by 
their momenta. We use the notation inrk.) = ik 1, k2, • • • k„) for the Fock states of both the Rindler and Minkowski 
fields. Note that we leave off subscripts denoting the Rindler and Minkowski Fock states and let the field operators 
imply the label, i.e. j'l'j)^ —>• |4'j) and jk^)^ —>• jk^). Working in the interaction picture, the acceleration-induced 
probability amplitude for our massive initial state to transition into JV = tir + n final state particles is given by. 


n nji 

e=i j=i 

The magnitude squared of the transition amplitude gives the differential probability per unit momenta of each 
Minkowski particle, i.e. = |Mp. Note that we are using the more compact notation ]Xj=i Moreover, 

since we are taking the magnitude squared of a complex integral, we remind the reader that there are two dummy 
variables in expressions such as [//(x)dxp = ff dxdx'\f{x)\‘^, where |/(x)p = f{x)f*{x'). Thus the differential 
probability of our V-particle acceleration-induced transition is given by 
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( 4 ) 


By factoring out the n complete sets of momentum eigenstates we can further simplify the above expression. The 
resultant two-point functions, i.e. Wightman functions, characterize the probability for each Minkowski particle to 
propagate along the accelerated trajectory. The product of the n Wightman functions then characterizes the total 
probability that all of the n Minkowski particles simultaneously propagate together. Moreover, we will endow each 
Wightman function with the index m to label each of the Minkowski particles and also the relevant observables 
computed. Hence 


V = 
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( 5 ) 


The inner products over our Rindler fields in the above expression select the appropriate mode function of each 
particle. For Rindler particles [12], these mode functions can be written as /^[a;(T)]e““'^. Analyzing the transition 
process in a frame comoving with the initial accelerated particle allows us to parametrize the system with its proper 
time. This implies that its energy is just its rest mass, i.e. uji = rrii, while the final particles have an arbitrary energy 
ujr- The inner products then imply 




"K 


j=l r=l 


hi [x{T)]e [x{T)]e 

r—1 

riR 

hi[x{T)] h^[x{T)] 


AEt 


( 6 ) 


Note that we have defined the Rindler space transition energy AE = — rrii to be the total energy difference 

between the final and initial Rindler states. We then find our acceleration-induced transition probability, Eq. (5), to 
be 


V = 
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( 7 ) 
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We can now use the overall normalization associated with the overlap of the spatial waveforms k [4,8] to simplify 
the above expression. In experimental settings the initial accelerated particle will be described by a wave packet 
with mode functions of the form />i>(a:) ^ The spatial distribution of the initial particle, in the 

directions perpendicular to the acceleration, is assumed to be finite, e.g. Gaussian, and is given by With the 

mode functions properly normalized [13], k will be of order unity. Hence 


= / / 


hi[x{T)] /ija:(r)] 


r—1 


( 8 ) 


Having integrated over the spatial coordinates of the right Rindler wedge and noting that we still have to integrate 
over the proper times, it should be mentioned that the trajectories which characterize the Wightman functions must 
only depend on the proper time and have no spatial dependence. With this in mind, our transition probability 
becomes 


V = G 


2 

n 


Y[ G^[x',x]. 

m—1 


(9) 


Upon inspection of the above transition probability we find that we have now effectively reproduced the formalism 
the would be obtained if we had used an Unruh-DeWitt detector provided we identify the energy gap of the detector 
with the initial and final state Rindler energies AE. To better incorporate the time dependence of our current analysis 
we will make the following change of variables to the rapidity variables u' and u defined by u{t) = a{f)df. Hence 


cItcIt' 


dudu 

dudu 


, dr dr' 

du du' 
1 1 


( 10 ) 


We dropped the proper time dependencies in terms of the more compact notation a' = a(r') which we also apply 
for all other variables. In terms of the rapidity, the transition probability then takes the following form 


'P = Gl j j drdr'e ^ Gt,[x',x] 

m—1 

= Gl 11 dudu’-^e-^^^G'-P fr Gi[x',: 

n n ' 


( 11 ) 


Finally we decouple the integrals by formulating the new variables in terms of the difference ^ = u' — u and average 
V — " 2 ^ rapidities. The inversions of these transformations are given by u' = ^/2 + ry and u = —^/2 + ry. With this 
new rapidity parametrization, the transition probability takes the similar form 


V = Gl 


dr^d^-^e-^^^G'- 

a' a 




m—1 


( 12 ) 


For the sake of clarity we note that under this parametrization the primed and unprimed variables will all have 
dependencies such as a' = a(^/2 + 7y) and a = a(—^/2 + 7y). All components of the integrand depend on these variables 
in this manner. Moreover, we will eventually want to determine the transition rate, i.e. the probability per unit 
time F = ^. We shall choose the proper time parametrization Tr^ that characterizes the rapidity variable ry via the 
definition dry = a{Tri)dTri for this purpose. Finally, using the notation a(T,,) = a,, we find the transition rate to be 


P = 


Gt 


drjd^——e 
a a 


— iAE{r' — r) 


n Gt[x\x] 

m—1 



n Gti[x',x]. 

m—1 


=> F 


(13) 
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We mention as well the ability to compute the differential transition probability per unit rapidity ^ that 

follows along with this derivation as well. To develop the integrand into a more useful form, we Taylor expand the 
proper time interval about the point ^ = 0. Recalling the form of the coordinate transformations used above, we find 


To 


dr' du' 
du' d^ 


^ - To + 


?=0 


dr du 
du d^ 


«=o 




Qi^i 2 


Qj'q 2 


(14) 


Similarly, we expand the proper accelerations about the same point but we also disregard terms of order jri/a^ and 
higher. Hence 


a '(^/2 + ry ) 


da' dr' du' 
dr' du' d^ 




T/ 1 1 






Qfj [1 

Qr,. 


a2 2 


(15) 


Similarly, we obtain a{—^/2 + ry) ~ a,; as well. We must also recall that all other components of the integrand must 
be expanded to the appropriate order. As such, we employ the notation Gy^[x',a;]r, to imply the necessary Taylor 
expansion of the Wightman function. Moreover, now that acceleration is a function of one variable only, we will drop 
the rj subscript and define = a(T^) = a. Thus our generalized response function takes the following form 


1 r ” 

r = Gl- n ( 16 ) 

^ J 1 

This expression determines the transition rate for a Rindler, i.e. accelerated, particle to decay into tir Rindler, 
i.e. accelerated, particles accompanied by the simultaneous emission of n Minkowski, i.e. inertial, particles. This 
is the process that is expressed schematically in Eq. (1). A conceptual depiction of the above acceleration-induced 
process is also illustrated in Fig. 1 for clarity. To finalize this section we comment on the relationship between the 
method of field operators, used in this manuscript, with the alternative method of detectors [8]. The method of fields 
enables a more general analysis via the inclusion of an arbitrary number of Rindler particles of varying energy in 
the final state. The energy difference, as measured in the proper frame of the initial accelerated particle, is given by 
/S.E = J2r=i‘^r — iTii', see e.g. Eq. (6) and the subsequent discussion. That is, the difference between the sum of all 
final state Rindler particle energies and the initial accelerated particle’s mass. To map the analysis to the method of 
detectors, all one needs to do is identify this energy difference with the energy gap of the detector, i.e. a two-level 
system. In short, the initial and final state energies of the Rindler particles, as measured in the proper frame of 
the initial accelerated particle, define the two energy levels of the detector. For a more in-depth discussion of this 
correspondence we refer the reader to Ref. [8]. In the next section we develop the time-dependent formalism that will 
be used to compute the Wightman functions and their subsequent Taylor expansion used in Eq. (16) above. 
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FIG. 1: A pictorial representation of the acceleration-induced transition from Eq. (16). 


III. TIME-DEPENDENT SPACETIME TRAJECTORIES 

Prior to evaluating the Wightman functions we will need the time-dependent spacetime trajectories the Minkowski 
fields propagate along. For the time-dependent proper acceleration a(r), we recall the rapidity is defined by u{a) = 
a{T)dT. Using the Rindler chart to characterize the resultant accelerated motion, see e.g. [14], we have 


x{t) = 

j da sinh [«(«)] + Xq 

t{T) = 

j da cosh [u(a)]-I-to- 


(17) 


In order to compute these, and related, integrals we shall make use of a simple variant of the method of steepest 
descent. Here, we will maintain the appropriate expansion and again disregard terms of order j/a^. As such, the 
method utilized can accommodate general acceleration profiles quite easily but depending on the system additional 
care must be taken if the acceleration goes to zero within the interval. Thus we consider the following change of 
variables and integration by parts [15], 


I = 
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r' 1 d 

T ±a(a) da 
1 




o±“(“) 
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o±“(“) 
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(18) 


Note that we used the fact that u'{t) = a(r) and u"{t) = j(r). Moreover we see that we obtained a solution to the 
integral to zeroth order in j ja} as required to be consistent with the development of the generalized response function 
in the previous section. We then have the general form of the integral to be applied to our spacetime intervals. Hence 




±a(T') 


^±u{t') _ 


o±“(t) 


±o(r) 


(19) 


We now apply the above formula to compute all necessary spacetime quantities associated with our generalized 
hyperbolic trajectory. Considering first the spacelike interval Acc = x' — x, we find 
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Ax = j da sinh [u{a)] + xo — da sinh [«(«)]+ 


da sinh [u(a)] 


1 r 


da[t 


Lt(a) _ -u{c 


n 


Y cosh - yt cosh [m(t)] . 

a{T) 


( 20 ) 


Utilizing the same coordinate transformation, i.e. u,u' —>■ ^,77 and Taylor expansion, i.e. a'(^, 77 ), a(^, 77 ) —>■ = a 

as in Sec. II, the above spacelike interval becomes 


Ax = 


-ttT cosh [u(t')] - -rw cosh [m(t)] 
a(r') o(t) 

cosh [5/2 + 77 ]- cosh [-5/2 + 77] 


a'(77,5) 




1 


= -(cosh [ 5/2 + 77] — cosh [— 5/2 + 77]) 
2 

= - sinh [ 5 / 2 ] sinh [ 77 ]. 


( 21 ) 


Note the use of the hyperbolic double angle formula to obtain the last line. Similarly, for the timelike interval 
At = t'-twe find 


At = 


da cosh [7t(a)] + to ~ dacosh[u(a)] 

■'t' 

da cosh [7t(a)] 

J da[e“(“) + 


^0 


-ttT sinh [u(t')] - sinh [^(t)]. 
a(r') a{T) 


The subsequent Taylor expansion to zeroth order in j/a^ then yields 


( 22 ) 


At = sinh [7((r')]-^ sinh [7t(T)] 

a(T') 

= -77-TT sinh [ 5/2 + 77] - 7 ^ sinh [-5/2 + 77] 

a'(v,0 a(?7,5) 

= -(sinh [5/2 + 77] — sinh [—5/2 + 77]) 

= - sinh [5/2] cosh [ 77 ]. (23) 

It should be noted that the addition of the complex regulator to the timelike interval can be added in without 
changing the computation. Next we compute the complex regulated square of the spacetime interval, Ax^ — (At —te)^, 
more commonly found in Wightman functions. We can evaluate this quantity by substitution of the previously 
computed spacelike and timelike intervals; however, in the interest of further developing the time-dependent formalism, 
we shall first carry out a few more manipulations to obtain a more general form. Thus 



— {At — ief' = [x{t') — a;(r)]^ — [t(r') — ^(t) — ief 



2 



/ dasinh [u(a)] 

— 

/ dacosh [u(a)] — 

/ da cosh [u(q;)] — ie 


1 ^ 

r' 

— 

/ da cosh [u(a)] — ie 


dasinh [w(a)] 

T 

p pt pt 

= // dad/? (sinh [«(«)] sinh [u(/3)] — cosh [«(«)] cosh [m(/ 3)]) + 2ze / dacosh[M(a)] 

= ~ JJ dad[3 cosh[u{a) — u{/3)] + 2i€At. (24) 

Note in the last line that we used the hyperbolic double angle formula and rewrote the term on the complex 
regulator as At. Focusing on the integral, we break up the hyperbolic argument into exponentials to obtain a more 
convenient form. Hence 


^u(a)-u{0) _|_ g-(ti(a)-ti(^)) 


JJ dadp cosh {u{a) — u{l3)) = ^ JJ dadjde^ 

= \J^ J^ df3e-<l^'> + ^J dae-“(“) J 

= / (25) 


Note in the last line that we interchanged indices a O /3 to combine the exponentials into one expression. We 
should also note a similar expression was obtained in [10]. Thus the square of the spacetime interval, along with its 
complex regulator, can be expressed in the following more compact form: 


Aa;2 - {At - ief = J J + 2ieAt. 

In turning our task to evaluating the above integrals, we find 
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(27) 


Finally, performing the necessary change of variables and Taylor expanding the acceleration we obtain the Hnal 
form of the ie regularized square of the spacetime interval. Hence 
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Note that, in the third to last line, we pulled the positive definite factor out of the e and in the last line 

we absorbed the magnitude of the acceleration into the regulator. We keep the overall sign of the acceleration so as to 
not change the direction of the shift along the imaginary axis by changing the sign of e. To finalize the development 
of the components necessary to compute each of our Wightman functions, we also require the use of a Lorentz gamma 
to boost forward and backward between the proper and inertial lab frames. The Lorentz gamma can be computed by 
taking the derivative of the inertial time. We also parametrize our proper time here using the variable Tri such that 
u{Tri) = rj as in the previous section. Recalling dt = jdr we have 


—t - — n 
d-Tn dTr, J 

7 = cosh [rj]. 


da cosh [zz(a)] + to 


(29) 


In the next section we will utilize the above dictionary of formulas to evaluate the Wightman functions and its 
variants. The Wightman functions will then be used to compute the acceleration-induced transition rate, power 
emitted, and spectra of emitted particles. 


IV. THE WIGHTMAN FUNCTION AND ITS VARIANTS 


In this section we compute the variants of the Wightman functions used for the computation of various observables. 
We will be working under the assumption that all Minkowski fields are massless. To evaluate the vacuum to vacuum 
two-point function [16], we use the canonically normalized field operator ^ = / (^ 2 Tr)^/W 2 ui +al.e~d^'^~^'^'>]. 

Thus, 
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(30) 


Since we are dealing with the emission of Minkowski particles by an accelerated field, it is advantageous to boost 
the momenta into the frame that is instantaneously at rest with the accelerated field. The integration measure is 
Lorentz invariant, i.e. When boosting the momenta forward into the accelerated fields’ instantaneous rest 

frame via k^ A/k'' = k^ , the Lorentz scalar nature of the exponent necessitates the boosting of the spacetime 
interval back via Ax^ —>■ {h~^)/;Ax'^ = Ax^. It should be noted that each Ax^ is a proper quantity that has just 
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been boosted to the same velocity but in the opposite direction and therefore all of the relativistic effects, e.g. length 
contraction, remain the same. Also, here and throughout, we will denote all proper quantities with a tilde. Therefore, 
our Wightman function, evaluated in the proper frame, is given by 


Gt\x' 


1 

2(27r)3 


Cj 


(31) 


Integrals of this form are best evaluated in spherical coordinates. Without loss of generality we rotate the coordinate 
system until the momentum lies along the z-axis and then we enforce the condition that our Minkowski fields are 
massless, i.e. ui = k. As such, the integral reduces to 
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The above Wightman function will be used to compute the acceleration-induced transition rate. To properly 
evaluate the integral we require an infinitesimal shift along the imaginary time axis to regulate the oscillations at 
infinity. This is accomplished by letting At —>■ At — ie where e > 0. Thus, 
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Note in the last line we used the Lorentz invariance of the scalar product to boost the spacetime interval forward 
into the lab frame. This will facilitate later computations and also highlights the appropriate Lorentz invariance of the 
Wightman functions. It serves to also remember the presence of the complex regulator within the time component of 
the interval. Using the appropriately Taylor expanded spacetime interval derived in Eq. (28) of the previous section, 
we find the Taylor expanded Wightman function to be 


G^[x',x]r 


I I 

(27r)^ Ax^^Ax^ 

I 1 

(27r)2 ^ sinh^ [^/2 - sgn {a)ie\ 
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(34) 


The Wightman function characterizes the probability for a particle to propagate along a given trajectory and is 
summed over all energies, i.e. integrated. Thus, by multiplying each probability by the energy (see e.g. [6]), and then 
integrating, we can effectively compute the average energy carried by the particle during the transition process, i.e. the 
power radiated. Denoting this energy weighted Wightman function as Qj^[x',x\r^, we then carry out its computation 
in a similar manner. Beginning with Eq. (32), we find 
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1 At — ie 

2i7r2 [Ai^ — {At — ie)^]^ 
1 At — ie 
2i7r2 [Ai^Ax^]^ 

1 Atj^ — ie 
2i7r2 [Ax^Ax^]^ 


(35) 


Note that the presence of the proper frame timelike interval in the above quantity necessitated the use of the 
Lorentz gamma to boost it back to the lab frame via At = At = At/q. Evaluation of this energy weighted Wightman 
function along the Taylor expanded hyperbolic trajectory yields 


1 At — ie 
2i7r2 [Ax'^Ax^]^ 

1 I sinh [^/2] — ze 

^ " 2 i ^2 sinh2 [^/2 - sgn {a)ie]]^ 
sinh [^/2] — z|e 
z(47r)2 sinh'^ [^/2 - sgn (a)ze] 
aP' sinh[^/2] — z sgn (a)e cosh [^/2] 
z(47r)2 sinh^ [^/2 - sgn (a)ze] 


z(47r)2 sinh^ [^/2 - sgn (a)ze] 


(36) 


Again we pulled a positive definite factor 2 cosh [^/2] out of the e. We also absorbed the magnitude of the acceleration 
into the regulator and then recombined the numerator into the hyperbolic sine since it was the Taylor expansion about 
small e. Finally, examining how the propagation, and thus emission probability, changes as we vary the energy enables 
us to compute the spectra. Taking the derivative of the Wightman function, Eq. (32), yields 


Ar± 

duj 


1 


1 


X , X = 




2(27r)^ iAx I 

Substitution of the relevant trajectories, along with the appropriate expansion, yields 


(37) 



auj 


1 


1 


2(27r)^ iAx l 

^ — iCbAt 


^iijj{Ax — At) ^—iCj{Ax-\-At) 


1 


(27r)2 Ax 

1 g-i‘iAi/7 


sin (wAx) 
sin (wqAx) 


UJ 


(27r)^ qAx 

a g-*“| sinh (5/2) 

(27r)2 sinh (5/2) sinh {2ri) " a 


sin I — sinh (5/2) sinh {2r]) 


(27r) 


(38) 


In the last line we expanded the arguments about small 5- Note that we have kept, to first order, a dependence in 
the phase so we still encode the dynamics. Now that we have our catalog of the Wightman function and its variants, 
the remaining sections will be devoted to the computation of the relevant observables of the theory. 
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V. TRANSITION RATE 

The generalized response function developed previously will now be used to compute the acceleration-induced 
transition rate. This will enable us to analyze the decay of unstable particles as well as the excitation of stable 
particles into states of higher energy. We begin by recalling the functional form of the response function, Eq. (16), 


1 f ” 

r = Gl- Gi[x',x]^. (39) 

m=l 

To compute the transition rate we use all n of the Wightman functions in the above product over our final state 
particles. Recalling the form of the Wightman functions, Eq. (34), we find that the transition rate of n-particle 
multiplicity simplifies to 


r = 


Gi- / d^e 

./ 


-iAE^/a 


n 


5 


m—1 


-iAE^/a 




( 47 r )2 sinh^ [C/2 — sgn {a)ie] 




^ — iAE^/a 


sinh^^ [^/2 — sgn {a)ie] 


(40) 


At this point we must focus the integral. Specifically we must note how the integration depends on the sign of the 
acceleration. By letting C —>■ sgn (a)C we see the integration is independent of the sign of the acceleration. As such 
we set sgn (a) = 1, leaving the transition rate dependent only on the magnitude of the acceleration |a|. Hence, 


r 



^-zAEi/lal 

^sinh^" [C/2 — ie] 


(41) 


We now see that the removal of the ie yields a singularity structure with poles of order 2n at C = 27ria with the 
integer cr > 0. Now that we know the pole structure, we rid ourselves of the complex regulator. Rewriting the 
denominator of the integrand in exponential form will yield a more useful form. Hence, 


r 


Gl 


filV" 1. /dr 

\47rJ jaj 7 sinh^" [C/2 — A] 

/mV” 1 f 

y [e«/2 - e-«/2]2» • 


Employing the more convenient change of variables w = we find 


(42) 


r 




27r J 


G 


2 ( 


la 
V 27r 


2n 


dw 


^-tAEi/\a\ 

[e?/2 _ e-«/2]2" 

^—iAE^/\a\+n—l 


[w- 1] 


2n 


(43) 


— •i/3 + 'Y 

The above expression can be integrated using the residue theorem. The integral is of the form with the 

appropriate identification for j3, 7 , and 5. Also note the added conditions that 7 is an integer and 5 is an even integer. 
Evaluation of this generalized integral yields 
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/ 


dw 


w-iP+'i 
[w — 1 ]*^ 


2TTi ^ 

(5 — 1)! dw^~^ 

2m r(l-i/3 + 7) r ^+^_5+i1 

(5-l)!r(-i/3 + 7-(5 + 2) ^ L 

27rz r(l-i/3+7) 27 i-cr/ 3 +i 27 r(T( 7 - 5 ) 

((5-l)!r(-i/3 + 7-(5 + 2) 

2m r(l — i /3 + 7 ) 1 

{6 - 1)! r(-i/3 + 7 - (5 + 2 ) 1 - 
27ri r(i/3 — 7 + (5 — 1 ) 1 

(5-1)! r(i/3 - 7 ) - 1 ■ 


[w- 1] 


5 W 


— i/3+7 


[zzi-ir 


(44) 


Note that we used the identity r( 2 )r(l — z) = 7 r/sin( 7 rz), along with the properties of 7 and S, in the last line. 
The sum in the above equation converges for < 1 which is true for negative /3. This corresponds to AE < 0, 
i.e. decays. In order to evaluate the sum we may assume /3 to be negative to yield the closed form expression of the 
convergent sum. We can then relax the condition for AE > 0 and we note that, in the zero acceleration limit, the 
transition rate appropriately diverges. This merely illustrates the fact that inertially stable particles have an infinite 
lifetime. Utilizing the above derived formula, we evaluate the integral in Eq. (43) to be 


/ 


dw 


yj-iAEi/\a\+n-l 

[w — 1 ] 2 ” 


The transition rate is then given by 


2m r(iAU/|a|+n) 1 

(2n - 1)! r(iAU/|a| + 1 - n) - 1 ' 


(45) 


p/Ap 27ri r(iAU/|a|+n) _ 1 

1 a{2n-l)lT{iAE/\a\ + l-n)e^^^^/\<^\-1 


(46) 


We will take one final step in simplifying the above gamma functions using the Pochhammer symbol for rising 
factorials (a;)“ = = n^rQ(a; + j). This serves to better present the resultant polynomials of multiplicity. 

Applying the Pochhammer identity to the above expression, along with the identifications x = iAE/\a\ + 1 — n and 
a = 2 n — 1 , yields 


r(fAU/|a| + n) 
r(zAU/|a| + 1 — n) 


2n—2 

(iAU/|a| + l-n+j) 

i=o 

n—1 

{iAE/\a\+k) 

k— — {n—l) 


n—1 


iAE 


J^(iAU/|a| + k){iAE/\a\ — k) 


k=0 


a 


n—1 


i-^r^Uii^E/ar+k 

k^O 

2n—1n—1 - 

n 

fc=0 L 


iAE 
fiAEV 

1 kl J 


i^r 


(47) 


Utilizing the above expression, along with the double factorial identity (2a;)!! = 2^(a;)!, we finalize the computation 
of the acceleration-induced transition rate. Thus 


r{AE,r^) = Gl 



1 

(4n- 2)!! 


n—1 


n 



1 

g2'7rAE^/|a| _ ^ ' 


(48) 
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The above transition rates are characterized by a bosonic distribution along with the integer indexed polynomial of 
multiplicity. Normalizing the rates to the multiplicity-dependent coupling via T = T/G^, we explicitly compute the 
normalized rates for the first five multiplicities. Hence 


ri(AG,a) = 
f2(AG,a) = 
f3(AG,a) = 
f4(AG,a) = 
f5(A£;,a) = 

The n = 1 case is the standard result for computing the transition rate of an Unruh-DeWitt detector. We should 
also mention that the acceleration-dependent lifetime is easily computed via r = l/T. We plot (see Figs. 2 and 3) the 
transition rates for acceleration-induced decay and excitation respectively. The fact that the excitation rate rapidly 
goes to zero, in the inertial limit, reflects the infinite lifetime for such processes. Moreover, it should be noted that 
there is an acceleration-dependent crossover scale where one multiplicity dominates the transition process. The next 
section deals specifically with these crossovers. 


AE 


1 


27r - 1 

487r^ g27rA£^/|a| — ^ 


AE^ l + 5{^)^+4{^y 


38407r5 - 1 

1 + 14(^)V49(^)S36(^)^ 

6451207r7 g2,.Ais/|a| _ i 

AE^ f + 30 + 273 + 820 + 576 


1857945607r9 


g27rA£^/|a| _ ^ 


(49) 



FIG. 2: The normalized transition rates, Eq. (49), with o = |a|/Ai? and AE = —1. 
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FIG. 3: The normalized transition rates, Eq. (49), with a = \a\//S.E and AE = 1. 


VI. MULTIPLICITY 

In this section we compute the dominant multiplicity as a function of acceleration. Specifically, for final states which 
contain either n or m Minkowski particles, we ask when the transition rate for an n-particle final state is greater than 
an m-particle final state. This will enable us to characterize the acceleration scale which selects a specific transition 
when there are multiple decay pathways. Let us begin by computing the inertial limit of the decay rates. This enables 
us to fix the couplings to that of the inertial decay rate A„. For decays we have /S.E < 0 and taking the limit 
a —>■ 0 of the transition rate, Eq. (48), yields 


An 

GS 


= G 


2 


1 


\ TT J (An — 2)!! 

/ 7 ^ \2n-l 

A„(^) (ta-2)!l. 


(50) 


In terms of the dimensionless acceleration d = |a|/AE, the crossover scale, i.e. when r„ = Tm, can then be 
computed in terms of the inertial rates. Hence 


An n [1+ n 

fc=0 j=0 


+ fa^]. 


(51) 


The ratio of the inertial decay rates Xm/^n is equivalent to the ratio of the branching fractions Brm/Brn = g oi 
each decay pathway [17]. Assuming n > m we find 
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n—1 

A„ I]^ [l + 

fc=0 


n [1+ kV] 


m—1 


Am n [l + 

i=0 
m— 1 


all [i+i"»“l 


n—1 

[1 + fc^a^] = g. (52) 

k—m 

The above equation defines the acceleration scale at which an n-particle multiplicity final state will dominate the 
TO-particle multiplicity final state. We now examine in more detail the cases when dk ^ 1, dk 1, n = m + 1, and 
n = m + 2. For the case of large acceleration, i.e. dk ^ 1, we find 


9 

9 

9 




n—1 




n [ 1 +^ 

n— 

n‘ 


k—m 

n—1 

~2(n-m) TT ;,2 


k—n 


-2(71—m) 


{m — 1)! 


(n — 1)! 

(m — 1)! J 

2 l l/(2(n-m)) 


L 1)! J 


(53) 


In the case of small acceleration, i.e. afc <C 1, we make use of the properties of logarithms to simplify the compu¬ 
tation. Hence 


9 


In (ff) 


n [1+ kV] 

k—m 

/ n—1 

In n [1 + 


n—1 

ln(g) = ^ \n[l + k‘^d^] 

k—m 

n—1 

In (g) = 

k—m 


n—1 

g = 1 -I- ^ k^ 

k—m 


(54) 


The last line followed from Taylor expanding the resultant exponential. The sum of squares evaluates to X]fc=m = 
A(n — m) — i(n^ — m'^) + A(n^ — m^). Thus the acceleration scale for multiplicity transitions at low acceleration is 
given by 


ai 


9- 1 


r(n — to) — — ’m?) + \(n^ — rn^) 


(55) 


It should be noted that the above formula is most applicable for two decay pathways with nearly identical branching 
fractions, i.e. g ~ 1. The case when n — to = I can be used to characterize the emission of an extra photon in the 
decay process. We can trivially solve for the acceleration in this case. Hence 
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Ol 


' g- 1 
(n - 1)2 


(56) 


The case when n — m = 2 can characterize the emission of an additional particle-antiparticle pair during the decay 
process. This case can also be solved exactly. Hence 


g 

g 




k—n—2 

[l -I- (n - l)^d^] [l -I- (n - 2)^d^] 

A/[(n - 1)2 + {n- 2 ) 2]2 + 4(^ - l)(n - l)‘^{n - 2)2 - [(n 

2(n- l)2(n- 2)2 


1)^ -I- (n 



(57) 


The above acceleration scales can be used to fine-tune a system to select a preferred decay pathway and also, if the 
acceleration of the system is known, predict the relevant branching fractions of the system under study. In the next 
section we shall analyze the power radiated away by an accelerated particle. 


VII. POWER RADIATED 

In order to compute the power emitted by the fth particle Si we use the energy weighted Wightman function 
for that particle in the generalized transition rate from Eq. (16). Thus 


V = Gi- die 


— iAE^{a 


n 


m—l 


n—1 


S^ = ^n~ j 


m^i 


(58) 


Note that we separated out the energy weighted Wightman function of the ith particle. Recalling the explicit forms 
of the appropriate Wightman functions, Eqs. (34) and (36), the power radiated away by the fth particle simplifies to 


= Gl- f 

die-^AEi/a^ 

a J 


= Gl- f 

die-^‘^^Ga 

a J 


/ . \ 2n+l 

= -GIAtt 


n 

\Att I 


n—1 

m^i 


a? 1 

1 

z(47r)2 sinh^ [^/2 - sgn (a)ze]. 

(47r)2 sinh^ [C/2 — sgn {a)ie]_ 


-I n—1 


- / di- 


^—iAE^/a 


sinh^"'*'^ [i/2 — sgn {a)ie] 


(59) 


Again, with the same rescaling i —>■ sgn (a)^, the integral is invariant under the change of sign of the acceleration. 
Thus we let sgn (a) = 1 as before. Hence 


S^ = -GIAt: 


47r ) 


2n+l 


di- 


„-iAiSC/|o| 


sinh^"’’’^ [i/2 — ie] 


(60) 


Note that we find a pole structure similar to the transition rate in the absence of the regulator, this time with poles 
of order 2n -I- 1 at ^ = 2'!Tia with integer cr > 0. Using the same computational machinery as before, we employ the 
change of variables as the previous section, w = e^. Thus 
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2n+l 


1 


i|o|\ 

Itt J \a 

2n+l 


5. = -G^47r - de 


„-*AE4/|o| 


= -G^47r ^ rr MC 


1 


= -G247r 


z|a 

2tt j \a 

• I I \ 2n+l 

27r y 


1 


dw 


sinh 2 "+i [i/2-ie] 

g-iAE4/|o| 

[eC/2 _ g-C/2]2n+l 
— iAS/|a|+n—1/2 




[w - 1]2»+1 

Examining the above integral in a similar manner as in the last section we note that it is still of the form 


(61) 

£2!£1I 

[td — 


but with odd integer 5 and with 7 being an odd integer multiple of 1/2. Let us also note that 7 — 5 will also be an 
odd integer multiple of 1/2. The effect of this will be an alternating sign in the sum over the residues leading to a 
different thermal distribution. Hence 


dw 


yj-iP+l 

[w — 1]*^ 


27ri 


■E 


J(5—1 


[u>- 1] 


5 W 


— i/3+7 


[zc-ir 


(5 — 1 )! dw^~^ 

27ri r(l — i /3 + 7 ) 

((3 — 1 )! r(—i/3 + 7 — 5 + 2 ) 

27ri r(l — i /3 + 7 ) 

(5 — 1 )! r(—i/3 + 7 — 5 + 2 ) 

2 m r(l — i /3 + 7 ) 

(5 - 1 )! r(-i/3 + 7-5 + 2 ) + 1 

27ri r(i/3 —7 + 5— 1 ) 1 


E 

(T —0 
00 

E< 

<7—0 


— 2^+7—<5+11 


W 


^2'K7/3-\-i2'77’(7{y — 6) 


1 


(5—1)! r(i/3 — 7 ) +1 ^ ^ 

Again we used the properties of 5 and 7 to manipulate the gamma functions in the last line. Then, utilizing the 
above formula, the integral in Eq. (61) yields 


dw 


w 


-*A-E/|a|+n-l/2 27ri r(iAE/ |a| + 71 + 1/2) 


1 


+ - l]2"+i (2n)! r(iAE/|a| + 1/2 - n) + 1 

The expression for the power radiated by the ith particle is then given by 


S, = -GfdTT +d 




2n+l 


1 27ri r(iAE/|a|+77+1/2) 


1 


(63) 


(64) 


2tt' ) |a| (2n)! r(iAE/|a| + 1 / 2 - 77 ) + 1 ' 

We again use the Pochhammer identity = 0^=0 + /) the identifications x = iAE/\a\ + 1/2 — 77 and 

a = 2n. As such, the above gamma functions yield a cleaner expression. Thus 


r(fAE/|a| +77+ 1/2) 
r(fAE/|a| + 1/2-77) 


2n-l 

{iAE/\a\ + 1/2 - 77 + j) 

3=0 

2n-l 

n (7AE;/|a|+72) 

iodd^-i‘2n-l) 


2n-l 

(zAE/lal+72)(fAE;/|a|-+2) 

^odd — 1 



( 65 ) 
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Using the same double factorial identity (2a;)!! = 2^(a;)!, we find that the power radiated by the ith particle is given 
by 


5. = Gl4. (^) ^ 


n 

k=0 


1 + 


2fc + l 


f—V 

\aeJ 


1 


^2'KAE/\a\ _|_ 2 ' 


( 66 ) 


We see that the above power radiated is characterized by a fermionic distribution as well as the half-integer indexed 
polynomial of multiplicity. The root cause of this change in statistics is the fact that we had poles of odd integer 
order rather than even integer order as in the transition rate. The first few normalized power functions, S = S/G'^, 
are then given by 


Si{AE,a) 


52 (AU, a) 


SsiAE, a) 


Si^AE, a) 


27T e27rAS/|a| _|_ ^ 


9671^ g27rA£;/|a| _|_ 


A 1 I 35 / g I 259 / g I 225 / a 

1 - 1-4 tAgj 16 VAgj 64 \ AE) 

115207r5 ^244AE/\a\ I 

AE^ 1 + 21 + 

25804807r7 g2,rAi!;/|a| i 


11025 1 g 
256 \Ae) 


(67) 


We note that the n = 1 case, along with AE = 0 as determined in [18], reproduces the known a? dependence for 
the power emitted by bremsstrahlung. The plots of the power radiated for both transitions up and down in Rindler 
energy can be found in Figs. 4 and 5 respectively. Note that the power radiated diverges in the case of a positive 
energy transitions. This reflects that inertially stable particles do not transition up in energy and radiate energy away. 



FIG. 4: The normalized power radiated, Eq. (67), with a = ajAE and AE — —1. 
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FIG. 5: The normalized power radiated, Eq. (67), with a = a/AE and AE — 1. 


VIII. ENERGY SPECTRA 

Utilizing the same prescription, we now endeavor to compute the energy spectrum of the ith Minkowski particle 
emitted in the transition process. The transition rate from Eq. (16) characterizes the probability per unit time that 
a transition will occur. The differential transition per unit energy then characterizes how the probability of emission 
changes with energy, i.e. the spectra. Therefore we begin by taking the derivative with respect to the ith particle’s 
proper energy. Hence 


dr 

dui 


diS, 


d^e 


— iA.E^/a 




m—1 


= Gl^ [ f[ Gt,[x' 


m^i 


( 68 ) 


Again, we have separated out the relevant variant of the Wightman function for the observable we are calculating. 
Referring to Eqs. (34) and (38) for the Wightman function and its derivative respectively, the above energy spectra 
reduce to a more convenient form. Hence 


dr 

duji 


Gl\I f[ G^[x\ 

m^i 


Gi 


Gi 


1 




UJ 


(27r)^ ) 




a 


-I n—1 


de 


(47r)2 sinh^ [^/2 — sgn {a)ie] 

p-d^E+Cj)^/\a\ 


sinh 


2(n-l) 


K/2 - ie] 


(69) 


Note in the last line that we enforced the invariance of the above integral under a change in sign of the acceleration. 
We have encountered similar integrals for both the transition rate and the power radiated. It is worth noting 
that the difference is that the frequency variable which we are Fourier transforming with respect to has shifted 
via AE —)> AE + Cj. We can evaluate the n = 1 case quite easily at this point. Hence 
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duj 


n—1 


r 2 ^ ^ f - -*(AiS+cD)C/|a| 

G?^(5(A£; + w). 


(70) 


The above expression is merely a statement analogous to Fermi’s golden rule. Note that the presence of the delta 
function serves to enforce conservation of energy in the case of one particle emission. This implies that when there is 
only one Minkowski particle emitted that the radiated particle, as measured in an inertial frame instantaneously at 
rest with the accelerated particle, carries away the total change in the Rindler space energy. For the higher multiplicity 
cases we rid ourselves of the regulator and note the similar pole structure of order 2 (n — 1 ) when ^ = 27rzcr with integer 
CT > 1. Finally, we make the same change of variables w = to obtain 


dCj 
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{2tt) 

UJ 
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ia \ 
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2(n-l) 


(27r)2 \ 2 tt ) |a| , 






(27r)2 \2ttJ 
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Q-i(^E+Cd)^/\a\ 

sinh^^”-^) [C /2 - ie] 
2(«-l) 


[e «/2 _ e-«/ 2 ] 


dw 


w 


— i{AE-\-Cj) /|a| +n—2 


[w- 1] 


2(n-l) 


(71) 


The above integral can be evaluated using the integration formula from Eq. (44) provided we make the relevant 
identification of the indices /3, 7 , and S. Moreover, since the 7 is an integer and S is an even integer, the Pochhammer 
identity holds as well. Thus, by making the identification of n —>■ n — 1 from the transition rate, we may merely quote 
the final form of the spectra. Hence 
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Glib f AE + bb 
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2n—3 
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(4n — 6 )!! 
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1 


g27r(A£;+(Xj)/|a| _ ' 


(72) 


Here we find a bosonic thermal distribution with a chemical potential and a polynomial of multiplicity characterized 
by an integer index. It is interesting to note that the total change in Rindler space energy is identified as the chemical 
potential of the thermal bath. Recalling the above spectra gives the probability of emission per unit energy per unit 
time, we note that it needs to be normalized via = Af- This serves to scale the overall spectra and remove 
the effective differential time averaging. Below we write out the first few spectra normalized to the coupling via 
7 ^^ = A/” as in the previous sections. Hence 
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7^2 (AE, a, w) 

7^3 (AE, a, w) 
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(73) 


In Figs. 6 and 7 we show the spectra normalized by the transition rate, i.e. TV), of the zth particle emitted for 
both transitions down and up in Rindler energy respectively. It is interesting to note that regardless of the sign of 
the transition energy, we find Planck-like spectra for all multiplicities greater than 1. 
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FIG. 6: The normalized spectra, A/i = , with a = 1 and AE = —1. 



FIG. 7: The normalized spectra, A/i = 


i4^, with a = 1 and AE = 1. 

r aujj ’ 


IX. DISPLACEMENT LAW 


To better characterize the spectra, we can also look at the peak energy of the emitted Minkowski particle for each 
multiplicity greater than I. In the interest of determining the maximum of each peak via setting the derivative equal 
to zero, we can drop all prefactors and focus only on the energy dependence. Labeling the spectra polynomial of 


multiplicity Ad„ = Y\l=o 


I + fc2 ( 


AE-\-Cb 


we have 
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g27r(A£;+(i)/|a,| _ 

Then, by taking the derivative with respect to Cj and setting it equal to zero, we find 
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(75) 


Note in the last line that we defined the dimensionless parameter x = ^(Ai? + CC). Now we must evaluate the 
logarithmic derivative of the polynomial of multiplicity. Hence 
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(76) 


Then, combining the logarithmic derivative with the Eq. (75), we obtain our numerically solvable displacement law 
which allows us to determine the peak energy of the emitted Minkowski particles. Hence 
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e® — 1 
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(77) 


Then, in terms of the numerically solved x, we find that the peak energy is given by 


w = xE - AS. (78) 

27r 

Thus we have shown that the emitted particle’s energy, in the limit of high acceleration or zero change in the 
Rindler space energy, is directly proportional to the accelerated temperature, i.e. Cj = xta- This is in agreement 
with Wien’s displacement law but now with a more general transcendental equation to determine the displacement 
constant. Moreover, this establishes the fact that acceleration has an energy associated with it that is quantum 
mechanical in nature. By reinserting all relevant physical constants, we Hnd the energy of acceleration S^ is given by 


S 


a 


xah 

27rc 


(79) 


The acceleration dependence of the energy of the emitted particles would have a clear signature at sufficiently high 
accelerations. Advanced experimental systems could be coming online in the coming years that may be able to verify 
these effects [13]. 
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X. CONCLUSIONS 

In this paper we have established a computational framework capable of computing various observables of 
acceleration-induced particle physics processes. To better analyze physically realistic settings, a time-dependent 
formalism was developed to compute the spacetime quantities that go into the Wightman functions and its variants. 
These were then used to compute the transition rate, multiplicity, power radiated, energy spectra, and displacement 
law for accelerated decays and excitations of arbitrary final state multiplicity. We found the transition rate, power, 
and spectra are characterized by integer and half-integer indexed polynomials and thermal distributions of both 
bosonic and fermionic statistics. For the spectra, we found the total change in Rindler space energy plays the role of a 
chemical potential of the thermal bath. The displacement law for the spectra predicts the peak energy of the emitted 
Minkowski particles have a proper energy proportional to the accelerated temperature. This implies that there exists 
an energy, of quantum mechanical origin, associated with acceleration. 
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